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ABSTRACT 

Optical microcavities liave received mucli attention over the last decade from different research fields 
^D ranging from fundamental issues of cavity QED to specific applications such as microlasers and bio- 

J^ sensors. A major issue in the latter applications is the difhculty to obtain directional emission of light in 

^~; the far- field while keeping high energy densities inside the cavity (i.e. high quality factor). To improve our 

understanding of these systems, we have studied the annular cavity (a dielectric disk with a circular hole), 

where the distance cavity-hole centers d is used as a parameter to alter the properties of cavity resonances. 

We present results showing how one can affect the directionality of the far-field while preserving the 
<H uniformity (hence the quality factor) of the near-field simply by increasing the value of d. Interestingly, 

\^r^ the transition between a uniform near- and far-field to a uniform near- and directional far-field is rather 

abrupt. We can explain this behavior quite nicely with a simple model, supported by full numerical 
' ry-j' calculations, and we predict that the effect will also be found in a large class of eigenmodes of the cavity. 

O Keywords: Optical resonators, annular microcavities, microlasers, whispering-gallery mode, control of 

'■^^ directional emission, potential/dynamical tunneling 

O 1. INTRODUCTION 
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Optical microcavities [T] are known to allow the buildup of large energy density resonances (high quality, 
high Q-factor resonances) . This characteristic is in turn well suited for applications such as lasers (mirror- 
^j,,^ less microlasers [11 13]) and biological/chemical sensors (narrow resonance shift based microsensors [4, '5^). 

r^ Whereas regular shaped high quality cavities (disk, toroid, sphere) display uniform far-field emission, 

Oh many applications require both directional far-field emission and large Q value. It is therefore of central 

importance to design directional emission capable cavities that can also bear large quality resonances. 
,_^ This paper is the first of two contributions (see [Sj in these Proceedings). Here we will be primarily 

^ concerned with the parametrical control of the onset of non-uniform far-field emission from a uniform 

Cn Whispering Gallery Mode (WGM) inside an annular cavity. We will show that the loss of uniformity in 

^Q the far-field happens rapidly over a small range of the center-to-center distance d (see figure 1(a) ) while 

^"T the near-field behavior of the mode remains WGM-like. Since this phenomenon is strongly correlated 

with a decrease of the Q-factor of the unperturbed cavity, it suggests that another escape mechanism 
JiiJ must exist simultaneously with the usual potential barrier tunneling mechanism. We will then present a 

f~^ simple model showing the competition between two escape mechanisms: tunneling versus transition to 

,-H escaping channels. 
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2. GEOMETRY AND GENERAL CONSIDERATIONS 

We simplify the electromagnetic treatment of the passive cavity by assuming optically thin structures 
(almost 2D structures) embedded in a surrounding medium of refractive index Uq. This assumption 
allows for the separation of two polarization states (TM and TE). For simplicity, we will treat the TM 
case. The monochromatic electric scalar field ip (e~*"* time-dependence) satisfies Helmholtz equation 

VV(^,0) + w^(?',0)fcV(?'» =0 (1) 

where n is the refractive index of the electromagnetic system and k is the wavenumber in vacuum 
{k = 2tt/X). The numerical solution of this equation is carried through a scattering matrix formalism [7j 
and by the use of Smith's delay matrix [8|. The eigenvalues, tj, of the delay matrix may be plotted over 
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a wavenumber range constructing a delay spectrum as in figure 2(a) The resonances are seen as peaks 
of the delay spectrum. The delay tj of a given mode is related to its quality factor Qj by the relation 
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Solutions of equation M may be expanded for r > Rq, Rq being the maximal radial extent of the cavity, 



as 
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where Hm (•) are Hankel functions, {Am} the incoming wave coefficients and B^ = X^m' Smm'Am' the 
outgoing wave coefficients with S being the scattering matrix of the cavity. 

We may also treat this problem in a completely classical fashion, replacing the electromagnetic field 
by a free moving particle with intersections at the boundary where the total internal reflection (TIR) 
condition is met when 

\p\ = sinx > PTIR = rio/nc- (4) 

nc is the refractive index of the cavity and x is the incidence angle relative to the normal. Phase space 
is then constructed from the coordinates (pi — Si/Ro,pi — sin^i of successive impacts with the boundary 
(figures 1(b) and 1(c)). 
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Figure 1: Geometry and phase space of the annular cavity, (a) A circular inclusion of refractive 
index Uh and radius tq is placed a distance d from the geometrical center of the cavity; (b) Classical 
phase space coordinates, arc length s (alternatively angular position <j) — s/Rq) and p — sinx; ('^) Phase 
space constructed from the impact coordinates {(pi, Pi} for different initial conditions; it consists of 2 
disctinct parts, a regular region filled with whispering gallery trajectories and a non-regular region filled 
(mostly) with irregular trajectories. The 2 domains are separated by the geometrical condition \p\ — pnr = 
(rf + ro)/i?o- The classical emission (escape) region, \p\ < Ptik, is also indicated. 

Our choice of the annular cavity [S] is motivated on the one hand by the observation that the solution 
of (II]) for this structure exhibits both high Q WGMs and low Q/non-uniform emission modes. On the 
other hand, the classical phase space associated with this cavity proves to be well-separated between two 
completely regular regions and an almost entirely chaotic Non-Regular domain, |p| < pNR(rf) — {d+rg) / Rq 
(figure 1(c) I. More specifically, the chaotic region of phase space extends beyond the classical emission 



region of the cavity if pnr > Ptir and is completely included in the emission region if pnr < Ptir- We are 
then brought to think of a system with a mixed phase space incorporating characteristics of the perfect 
disk cavity (high Q, emission through potential barrier, pnr < Ptir) with those of a completely chaotic, 
e.g. the stadium cavity [10], (non-uniform/highly directional emission through transport, pnr > Ptir )■ 
Furthermore, the relative importance of the two characteristics appears to be controllable by the variation 
of pnrXc?) with respect to ptir- 



3. EFFECT OF PARAMETER d ON THE FAR-FIELD BEHAVIOR 

Wc set the physical parameters of the cavity as nc — 3.2 (typical semiconductor refractive index), 
rio — nfi — 1 {ptir = 0.3125), Rq = 1 and tq — 0.2Ro. We keep d as a control parameter. The delay 
spectrum of the unperturbed disk cavity shows a large resonance {ct ^ 8 x 10^i?o, c the speed of light 
in vacuum) for a WGM (11,1) at kRo ~ 4.499 (figure 2(a)). The notation {m,n) stands for (angular 
momentum number, number of radial nodes). Placing an inclusion at the center of the cavity (d — 0) 
has essentially no effect on the delay value and the resonance position. Increasing d, we observe that the 
resonance delay abruptly drops from its initial unperturbed value near d = 0.27i?o (figure [2(c) ) . 



In order to show that this observation is not the result of a large deformation of the near-field behavior 
of the WGM (11,1), we define a contrast measure 
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This quantity enables us to obtain the fraction of the |7tt,| components of the perturbed mode other than 
the initial pure mode JtooI at radial position r > Rq and inclusion distance d. It is understood as a 
measure of the non-uniformity of the field relative to the unperturbed WGM state. Two limits of the 
expression are of interest: the far-field (r — > oo) and the near field (r = _Ro)- For |mo| = 11, the results 
are shown in figure 2(c) 
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Figure 2: Numerical results for the dielectric annular cavity, (a) Delay spectrum for the disk 
cavity with a large (11, 1) resonant state at kRg « 4.499,' (b) A sharp decrease of delay value (Q factor) 
for mode (11, 1) happens near d — 0.27i?o/ (c) The contrast measure Q) expressing the non-uniformity of 
the field calculated in the far-field (r — > cx); plain curve) and near field (r = Rq; dashed curve). Note the 
rapid increase of Cn value in the far-field regime near d = 0.27i?o while the near-field behavior remains 
unaffected until d ^ 0.5Rq. 



We interpret the results of figure 2(c) as evidence that the potential-barrier tunneling escape mechanism 
remains unaffected by the deformation up to rf ^ 0.5i?o (the near field is roughly the same over the interval 
d = [0, 0.5i?o])- Furthermore, since the far-field undergoes a sharp change near d = 0.27i?o while the 
near-field remains unaffected by the inclusion, we may conclude that the non-uniform far-field pattern 
must be the result of a chaotic cavity behavior where transport in phase space is the dominant escape 
mechanism. The next section attempts to put these conclusions on sounder grounds. 

4. MULTIPLE SCATTER-AND-TRANSMIT/REFLECT MODEL 

A simple model [IT] has been designed to describe the field inside and outside the annular cavity. The 
total field is separated in intermediate local interaction fields (see [T^ for a similar description of the 



closed annular system). This model is similar with the multiple reflections treatment of the field in a 
Fabry-Perot resonator. The local fields inside the annular region d + ro < R'q < r < Rq are expanded 
as a superposition of incoming and outgoing waves just as ([3]) (see figure 3(a)[ ). An internal scattering 
matrix S' is then derived defining the interaction of waves with the internal effective boundary and the 
reflection and transmission coefficients characteristic of the interaction with the external boundary. The 
main difference to the Fabry-Perot resonator, apart from the geometry, is the conservative aspect of 
the interaction with the internal effective scattering matrix (S' is unitary) and, obviously, the mixing of 
angular channels that result from a non-diagonal S'. 

Our first observation is that, at resonance, for a high Q mode of principal component |?7Zo| in the annular 
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Figure 3: Multiple scatter-and-transmit /reflect model, (a) The effective internal scattering matrix 
S' relating the incident waves on boundary r = R'q to outgoing waves; (h) Close up view of a single round- 
trip. The two different escape routes out of the main components |mo| (transmission out of the cavity and 
transition to other angular momentum components) are indicated with dashed arrows; (c) Transmission 
(plain curve) and transition probabilities (dashed curve) for initial mode \niQ\ = 11. The probability of 
each mechanism becomes equal at d = 0.282i?o- 



region, iterations through the partial fields have little effect on the amplitude of the main component 
even after a large number of round-trip interactions. Of course this is inline with the high Q nature of 
the selected mode. This is true for the range of parameter d e [0, 0.5i?o] and indicates that the near-field 
is only slightly affected by the inclusion. If we now concentrate on one specific round trip (figure [3(b)[ ), 
we can evaluate the total electromagnetic power, Pmix oc X^imi^m 1^™ Pi ^"^ returning to the WGM's 
main angular channels JttioI following the interaction with the internal boundary at r 
electromagnetic power, Pout « ^ 
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in angular channels JttioI escaping the cavity through 
the external boundary aX r — Rq . These values, properly normalized to the incident power at the start 
of the sequence, define transition (mixing) and transmission probabilities. 

In fact, Pmix is the probability of mixing transitions to other low Q angular channels following an 
interaction with the internal boundary, while Pout is representative of the direct transmission hy tunneling. 
Both probabilities are plotted as a function of d in figure 3(c) It is rewarding to observe that the 
probabilities instersect at d ss 0.282Po, near the value where an abrupt transition is observed in the 



far- field Cn measure from figure 2(c) 



5. CONCLUSION 

This paper presents numerical evidences that a high Q WGM inside a dielectric annular cavity can 
be phase-space engineered to modify the far-field emission properties while retaining uniform near-field 
characteristics by simply controlling the position of a circular inclusion. The observed results indicate that 
a competition between two escape mechanisms takes place as the control parameter d is increased. These 
mechanisms have been identified as the transition probability to angular momentum channels other than 
the main component of the unperturbed WGM and the transmission probability through the external 
boundary of the cavity. While the latter remains constant over a large interval of d, the former rapidly 
increases and eventually overcomes the transmission probability. The value of d at which the probabilities 



become equal is in turn associated with a rapid decrease of the quahty factor and, accordingly, with the 
increase of non-uniformity of the far-field behavior. 

Having modified the far-field with a choice of parameter d, one may consider using the radius tq of the 
inclusion to optimize the directional properties of the emission. This further engineering is addressed in 
the companion contribution [B]. 
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